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Abstract 



Nuclear multifragmentation process can be viewed as a recombination of nucleons into 
D ■ clusters of various sizes. In a combinatorial analysis, various moments of cluster size dis- 



tribution appear to be quite simple in terms of canonical partition functions. This simple 
model can also describe a branching phenomena and a clusterization phenomena. The pos- 
sibility of applying this fragmentation model in describing a hadronization of quark-gluon 



£> ' system and a strangeness production is discussed. 

o 

g ! 1 Introduction 

A nuclear system at low excitation energy may be described as a Fermi liquid of nucleons 
which are strongly coupled each other through the mean field interaction. On the other hand, 
a hot nuclear system at a temperature higher than the nuclear binding energy per nucleon may 
be treated as a Fermi gas. Within the region where the excitation energy is near to the binding 
energy lie the possibility of seeing a liquid-gas phase transition in a nuclear matter. However, 
except in infinite nuclear matter, the number of nucleons in a nuclear system is too small to 



be considered as a thermodynamic system. Only up to few hundred nucleons are involved in 
a heavy ion collision. A finite nuclear system with the excitation energy near to the binding 
energy breaks up into many smaller pieces instead of forming an equilibrated system of either 
liquid or gas. This excitation energy region is called the multifragmentation region and the study 
of nuclear multifragmentation phenomena requires nonequilibrium considerations and studies of 
fluctuation and correlation due to the finiteness of the system. 

Based on a combinatorial analysis, a simple exactly solvable model describing this multi- 
fragmentation process has been developed Jl|, § . This simple model fits various data in nucleon- 
nucleus and nucleus- nucleus collisions || ffl. The fits are as good as the fits using more com- 
plicated fragmentation model based on the thermodynamic statistical model ||, Monte-Cairo 
calculation ||, and percolation model J/J. Considering a control parameter as a function of 
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temperature or time, this model can describe branching process [BJ and dynamical behavior of 
fragmentation or clusterization [[J. 

At much higher excitation energy, nuclear system may transform into a quark-gluon plasma. 
This transformation may be viewed as a fragmentation process of a nucleon into quarks and 
gluons. At low energy, quarks are strongly bound inside of a hadron while the asymptotic freedom 
for quarks may be applicable at a temperature higher than the hadronic binding energy. On the 
other hand, as the system cools down, quarks and gluons clusterize forming colorless hadrons. 
Both the formation of quark-gluon system and its hadronization have a similar structure with 
the nuclear fragmentation process. 

This talk summarizes a fragmentation model (Section |2|) and its branching and evolution 
aspects (Section |3|) which are discussed in more detail in Ref.@. In Section ||, the possibility 
of applying the simple fragmentation model in describing a quark-hadron phase transition is 
discussed by considering pions as a simple example. 



2 Fragmentation Model 

A method for investigating the fragmentation of a nucleus with A nucleons into pieces of 
smaller size is developed [p], 0, ||. The partitioning of A initial nucleons into groups of varying 
sizes gives rise to a distribution of cluster sizes. Let n& be the number of clusters of size k. The 
n\ is then the number of nucleons not in clusters. A partition 

™ = { n i} = {ni,n 2 ,n 3 , ■■ •) (1) 

defines a particular arrangement of the initial nucleons into clusters (112,113, ■ ■ ■) and monomers 
(ni) regardless of which nucleon belongs to which cluster. The number of clusters in a particular 
partition is called the multiplicity M; 

A 

M = Y^n l . (2) 

i=i 

A constraint exists which is A = Y^t=i n ih i-e., the conservation of the total number of nucleons 
in the nuclear multifragmentation process. 

Now, let us consider a more general case for the later use in Sections || and |[ Suppose we 
have a system with N different species, each species, say the species k, appearing times in 
the system. When each species k has an additive quantity the total value of this quantity of 
the whole system is given by 

N 

A = ^2rnai (3) 
i=i 

for a particular partition n. Here A may have any value depending on the «j. However, when we 
study the nuclear fragmentation phenomena in terms of the cluster sizes, «j = i is the number 
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of nucleons in each cluster and A is an integer representing the total number of nucleons in the 
system. We can also specify each cluster by their energy, i.e., a, = E iy and then A = E = J2i ^Ei 
becomes the total energy of the system. For q.{ = 1 as in a jet fragmentation process, a» and A 
are the number of jets and A = M. In the hadron-quark phase transition of nucleons, and A 
may represent the baryon number of quarks and nucleons. 

In a nonequilibrium phenomena, a system may have many different partitions n. For a 
collision involving 235 U -f 235 U, the number of possible partitions n is of the order 10 20 . Weighting 
each species i by parameters and we can assign a weight 



G A (n,x) 



T(A+1) 
Ml 



Ml 



N 



n , 



x 1 



N 



T(A+l)Yl 



i=i iri 



(4) 



ni\n2'. ■ ■ ■ n at! 

to each partition n. Here x = (xi,X2, £3, • • •) and the T(A + 1) is the gamma function which is 
Al for an integer A. The is the Boltzmann factor arising from the coexistence of identical 
species i in the partition. The parameter /3j may be related to the partition function factor of 
species i originating from its internal structure and Xi to the partition function factor related to 
the species i as a whole. In a nuclear fragmentation with = i, the choice of = i counts for 
the cyclic rearrangements of the nucleons in a cluster of size i, and Xi oc e~ Fl<yT ' v ^ kBT with the 
free energy Fi(T, V) of the cluster i at a temperature T and a freeze out volume V. 

Summing the weight Eq.([|) over all possible partitions with a fixed A, the weight for all 
the partitions having a fixed value A becomes 



iv 



Qa{x) 



53 G A (n,x) 



E r(A + i)n 

{Ma 
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Xa 



(5) 



This Qa{x) is then the canonical partition function of the system with a fixed value A. For the 
case of an integer A such as in a nuclear fragmentation, we can get easily the canonical partition 
function from a generating function 

M 



Q(u, x) 
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Ga (n, x) —— r 
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which is an infinite sum of multinomials. By summing the weight Ga{u, x) over the partitions 
with a fixed A and a fixed of species k, we obtain the weight 
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iV 



C^fn, x) 



E r(A + i)n 

{ni}A,n fc 
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X: 



ft** 



(7) 



for all the partitions having 77^ species of type k with a fixed value of A. 

The probability for a specific partition 77 to appear in the canonical ensemble (Qa) of a 
fixed value A is given by 



P A (n, x) 



G A (n,x) 
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We also define the inclusive multiplicity distribution 

P A (n k) x) = p A(n,x) 

{ni}A,n k 



Q n A k (x) 

Qa{x) 



(9) 



as the probability of having n k species of type k in the canonical ensemble. 

Due to the fit dependence of the weight G A (n, x) with the specific form of and (xi//3j) 
the derivatives of Eq.(f|) lead back to the weight GA{n,x) itself with a constant factor; 
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GA{n,x) = n^n^GAin-tX) 



(10) 
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G A (n,x) 
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(11) 



where n' is used to discriminate partitions having J2i n i a i — {A — {ma k + laj}) from the partitions 
n with A. Notice that Eq. ([T0|) becomes zero only when rij or n k is zero in contrast to Eq.(|ll|) 
which is non-zero only when I < rij and m < n k . Eq.([IT|) also applies to the case of j = k with 
rij = {n k — m), i.e., replacing the factorial factors on the right hand side by n k \/(n k — m — 
Due to the form of Eqs.(^) and (|TTD, it is quite simple to consider various average values in the 
canonical ensemble (Qa) of a fixed value A with the corresponding probability of Eq.(||): 
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Here ( ) represents an average over the canonical ensemble Qa{x) and 
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(13) 



is the m-th order moment of n k . Eq . (p!2|) also applies for j = k when rij = n k — m. For this case, 
Eq. fli~2"l) becomes Eq. ([T3l) with the m of Eq. (|l3"D set equal to m + I. The inclusive multiplicity 
distribution Pa(^/c, x) of a species k in a canonical ensemble can be looked at through a quantity 
defined by 

T(A + 1) 



Fa%i{x) 



(n k - m)\ 



Q n A k (x) 
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n k \ 
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for m < n k , which is related with the moments through Y A (k,x) = ELo^ra(^) / Qa(x). 
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Once we know the canonical partition function Qa{x) of a canonical ensemble with a fixed 
A, we can easily obtain various moments of the species distribution nk by Eg . fll2D . Specifically, 
the choice of aij = fa = i leads to various simple exactly solvable canonical systems || . When 
all the species have a same weight Xi = x, the canonical partition function simply becomes 
Qa(x) = T(x + A)/Y(x). This simple exactly solvable model describes the average behavior of 
nuclear multifragmentation phenomena very well H f|, [| [K|. The model were also able to 
describe fluctuations, correlations, intermittency, power law and scaling behavior involved in the 
nuclear multifragmentation phenomena |2|, |4], [11], |l^] . There are various sets of x each of which 
gives an exactly solvable model |J. A specific example among them is the case of Xi — (1 — q l )~ l 
which can describe the even-odd oscillating behavior appearing in various fields such as a pairing 
effect in a nucleus or the size of parties reserving a table in a restaurant . 



3 Evolution and Branching Phenomena 

A particular set of Xi values gives a specific distribution of possible partitions through 
Eq.(P) which in turn determines uniquely the various moments of the multiplicity n^. Now 
suppose the weight parameter Xj's are functions of an evolution parameter t such as time or 
temperature, 1.6. « % — j, 

(t). A change of the evolution parameter t of the system induces a 
change in the Xi(t)'s which assign a different weight given by G^in, x) of Eq.(|4]) for different 
partition n. This change in the weight occurring in response to the change of the evolution 
parameter t of the system gives the t-dependence of the partition function Qa(x) of Eq.(|]) which 
describes the evolution of the system in t, i.e., how the average distribution (n^) of species k in 
a canonical ensemble having a fixed value of A changes with t. 

The evolution of the canonical system can also be viewed as a redistribution of the quantity 
A = J2i n i a i m to various species k. The distribution of = n^a^ can be represented by 



N?(A) = a%Y?(k,2) = aT E F Am(x) / Qa(x). (15) 

n k =0 I 

(In Ref.0, N™(A) was defined by at}.Y™(k, x).) For the case of «j = i, = n^k is the mass 
distribution in a partition and N^(A) = (rrik) = k (nk) can be considered as the total average 
number of nucleons belonging to clusters of size k in the fragmentation of a nucleus of size 
A. The quantity N^(A)/A = {k/A) (n^) is considered in Ref.0 as a probability. In a nuclear 
multifragmentation process described in terms of the cluster sizes (a% = i), the evolution of a 
canonical system corresponds to the redistribution of nucleons into various size of clusters. On 
the other hand, when aj = Ei is the energy of a cluster i, the evolution is the redistribution of 



energy of the system into various clusters. Using Eq.([TT|) and 4 = Y.k=i 7h^ nx k(t) 
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Xk dx k 
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and separating j < k and j > k, we obtain the dynamic equation for N k l (A) as 

k-l N-k 

E J lk-i)( A ,k,rn) + E ^(fe+o^^m), 
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(17) 
(18) 



Oi(k+l) 

Eqs.(^6|) and (fL7|) show that the changing rate of N™(A) depends not on the mean behavior 
but on fluctuations and correlations in the distribution of m k = n k a k . The distribution is in 
equilibrium when Jj(A, k,m) = 0, i.e. when C l k = or when there are no fluctuations nor 
correlations. 

While Eq. ( |TED seems complex, this equation has a simple interpretation. The first term of 
Eq. (|T6|) represents the creation of species k arising from the fusion of one of the Nh* (A) species 
(k — I) and one of species / when the J l t k -n is positive (process (a) in Fig.|l] which is shown for a 
nuclear fragmentation process). If the Jh-i) ^ s negative, we can also interpret this same term as 
an annihilation process (b) of one of N^(A) species k. The second term of Eq.flT6|) represents the 
formation of species (k + l) arising from the fusion of species k and I (process (c) in Fig.|l|), which 
thus reduces species k, when J l k is negative. When J l k is positive, this second term represents the 
formation of species k arising from the fission of species (k + I) (process (d) in Fig.|]). 
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Figure 1: Evolution of fragmentation process in the fragmentation block diagram representation 
M. Each box represents a nucleon and the stack of boxes represents a cluster. 



For the m = 1 case, Eq.(^) represents the evolution of the average distribution (mk) 
®k ( n k) of species k 0; 
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where 



a 



in 



E 



(20) 



The result of Eq.fll9|) can be related to the Chapman-Kolmogorov equation || in a continuum 
limit of the species index j |fj]. Considering only / = 1 in Eq.(^), the result is further reduced 
to a Fokker-Planck equation || Due to the discrete value of k, Eq. (|l9"D is applicable for a 
finite system with the finite number of k values. 

The evolution of the inclusive multiplicity distribution P^n^, x) of a species k in a canonical 
ensemble can be looked at through a change of F^ m (x) defined by Eq. (|14D . The evolution of 
the F^ k m (x) contains information regarding the underlying dynamics such as in various cascade 



phenomena. Using Eq. ([□]), we obtain the dynamic equation for the function F^ m (x) of Eq.fJT 
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This result may be considered as a dynamical equation for the factorial distribution of species 
of type k in a canonical ensemble. The first term represents the dependence on the multiplicity 
n k of the same species k and the second term represents the dependence on the multiplicity rij 
of species of type j ^ k. The Eq. (^l|) for m — 1 is the changing rate of the number of species of 
type k, nk = F2\(x) / Q n ^{x), when each species i has the weight Xi(t). 



Initial : n^-jets 




> A -jets : Final 



Figure 2: Branching diagram ||. Each vertex (solid circle) represents a binary branching of 
Eq. (|2"3|) . The rik initial jets branch into the final A jets at (A — n/J-vertices. 



An interesting application of Eq.(|21|) is to branching processes and jet fragmentation. For 
— A = 1, in which Qa(x) = (Si=i x i) A — Y,n k ^lo(^) * s a simple multinomial @ with A = M, 
Eq.(|2T|) reduces to 
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Furthermore, if N = 2 with x\(t) = l/w(t) and x%(t) = 1 — l/u(t), then the Eq.(p2|) is reduced 
to a generalization of a dynamical equation (m = 1) for the multiplicity distribution function of 
jet fragmentation in 3 theory || having three lines at each vertex. For the species k = 1 with 
m — 1, the solution of Eq.(p2|) becomes the Furry distribution 



5» - fRr^TI) (zTi'-zT (23) 

in 3 theory with the evolution parameter uo{t) = e^. Here the A is the number of final jet 
lines at t and nk is the number of initial jet lines at t — 0. In 3 theory, each initial line (stem) 
branches into N = 2 new lines (branches) at each step (Fig.^). Here the species k — 1 is the old 
jet and the species j = 2 represents the new born jet from the old one with the probability of 
X2 — 1 — 1/uj. There are (A — n^) branching vertices. For a given A final jets, the initial jets 
should branch (A — n^) times producing (A — n&) new jets. is the probability of starting 
from nk initial jets to have final jets of a fixed number A. Detail discussions concerning the 
relation between the Furry distribution and branching processes can be found in Ref.||. 



4 Quark-Hadron System 

Study of quark-hadron phase transition is more complicated than the study of nuclear 
multifragmentation process. In a nuclear multifragmentation region, the nuclear system can be 
considered as a system constituted with one type of elements (nucleons) and the total number 
A of nucleons is conserved. In contrast to this, hadrons are constituted with quarks having 
different colors and flavors. The constituents carry color but the hadrons are colorless. We also 
need to consider anti-particles in this energy region thus having a non-conserved total number of 
particles. Further investigations are needed in order to use the fragmentation model to describe 
a hadron-quark phase transition. Here we just discuss how the fragmentation model may look 
like when we apply it to pions (simplest hadron). 




Figure 3: Fragmentation and hadronization of pion. 
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At high energy, pions may change to quark-gluon plasma phase. If we consider only quarks, 
then it would corresponds to a fragmentation of pion (tt) into a quark(g) and an antiquark (q) 
which is depicted by Fig.|3](a). A quark- antiquark system, in turn, clusterizes into a pion as shown 
by Fig.^|(b). This simplified case is quite similar to the nuclear multifragmentation phenomena 
discussed in previous sections. The number of particles are conserved for this case. 

To consider new particle creation, we need to consider gluons (G) too |14[], which may be 
created by the process of Fig.[3](c). Due to the self- interaction of gluon, the gluon may break up 
itself further giving more gluons, Fig||(d), which is the same form as the branching phenomena 
of Fig.§. After the system has cooled down finally, there are only pions without free quarks and 
gluons. Thus we also need the annihilation process of gluon (qq pair production), Fig.|3|(e). The 
processes in Fig.|3| may be viewed as a fragmentation process of Fig.|l] or as a branching process of 
Fig.|2|. The coefficient C l k defined by Eq.flilf) becomes and with Cg = C? for the process 
of tt -> qq, Fig.|(a). We also have C ? G = Cf and C q G = C%. Of course Q ^ C| here; the rate of 
process (a) is different from the rate of process (b) in general. For the process of Fig.Q(d), the 
branching parameter u(t) of Eq.(|23l) may be assigned. 

For the case of FigfJ, we have four different types of particles (species k); tt, q, q, and 
G. Initially we have n n (0) pions without any quarks or gluons, i.e., n q (0) = 0, n q (0) = 0, 
and nc(0) = 0. The hot pions break up into qq pairs according to the process (a), and then 
processes (c), (d), and (e) create more quarks and gluons. In the intermediate step, the system 
may have all four species, i.e., n n (t), n q (t), n q (t), and ric(t) all are nonzero in general with 
weighting parameters x^t), x q (t), x q (t), and xc(t) respectively. After the system has cooled 
down, we again have n q (oo) = n q (oo) = hg(oo) = with n n (oo) > n n (0) ^ 0. If we put 
x q (oo) = x q (oo) = xg(oo) = 0, then we are secured to have n q (oo) = n q (oo) = n G (oo) = 0. A 
more detailed study of this case would reveal a method of applying our simple fragmentation 
model to a quark-hadron phase transition including strangeness. 



5 Summary 

We have a simple exactly solvable model for nonequilibrium fragmentation phenomena 
based on combinatorial analysis. This model can describe nuclear multifragmentation and 
branching phenomena. The ubiquitous pairing phenomena is also able to be considered with 
this model. 

We have also looked at the similarities and differences of quark-hadron phase transition to 
our fragmentation model. But for a practical application of the model to quark-hadron system, 
we need further studies, such as what and we should consider, what Xk{t) should be, 
how many particles (N) we need to consider, etc. If we consider different colors and flavors, 
the number of species increases much more. Including the strange quark, we might be able to 
investigate the role of strangeness in a nuclear system. 
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